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Abstract  -  The  paper  presents  an  analysis  of  the 
phased  array  radar  allocation  demands,  when  track¬ 
ing  highly  manoeuvrable  anti- ship  missiles  (ASM) 
using  a  collocated  radar /IRST  sensor  combination. 
The  motion  of  the  ASM  is  modelled  using  the  quan¬ 
tized  acceleration  levels.  The  principal  aim  of  this 
analysis  is  to  determine  an  upper  bound  on  the  av¬ 
erage  radar  update  time.  This  bound  follows  from  a 
Cramer-Rao  type  error  bound  for  the  estimation  of 
linear  jump  Markov  dynamic  systems  [1].  Given  a 
dynamic  motion  model  of  an  ASM,  the  IRST /radar 
sensor  characteristics  and  a  tolerable  level  of  target 
state  estimation  error,  we  can  theoretically  predict 
the  maximum  average  update  time  required  for  the 
phased-array  radar.  The  presented  analysis  allows 
us  to  quantify  the  IRST  benefits  in  ASM  defence, 
without  a  need  for  extensive  Monte  Carlo  simula¬ 
tions. 

Keywords:  Tracking,  performance  bound,  resource  allo¬ 
cation,  anti-ship  missile  defence,  IRST,  pliased-array  radar, 
data  fusion. 

1  Introduction 

Modern  anti-ship  missiles  (e.g.  sea-skimmers,  BUNT 
missiles,  anti-radiation  missiles)  combine  low-altitude 
flight,  low  radar-cross-section,  high  speed  and  maneu¬ 
verability,  and  as  such  represent  a  serious  threat  be¬ 
yond  the  capabilities  of  a  radar-only  surveillance  sys¬ 
tem.  In  an  attempt  to  overcome  this  limitation,  a  mod¬ 
ern  shipboard  surveillance  system  would  typically  com¬ 
plement  a  phased  array  radar  with  a  passive  infra-red 
search  and  track  (IRST)  sensor  [2,  3].  Some  of  the 
characteristics  of  a  generic  IRST  sensor  are:  it  is  a 
passive  sensor  thus  it  cannot  in  principle  be  detected 
and,  consequently,  jammed;  it  provides  angular  mea¬ 
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surements  (azimuth  and  elevation)  of  target  position; 
it  is  not  affected  by  multipath  at  low  elevation  angles; 
it  has  a  very  high  angular  accuracy  and  resolution  wrt 
radar;  its  detection  range  against  anti-ship  missiles  is 
speed  dependent  (longer  range  for  higher  speed)  and 
consequently  provides  a  relatively  constant  “time  to 
closest  point  of  approach” . 

The  principle  of  operation  of  an  IRST  aided  radar 
surveillance  is  as  follows:  the  IRST  passively  scans  the 
horizon  at  a  constant  scanning  interval  in  order  to  de¬ 
tect  low  altitude  threats.  Each  such  detection  serves  as 
an  alert,  that  is  then  used  to  allocate  and  cue  an  agile 
beam  confirm  dwell  with  a  pulse-Doppler  waveform  of 
much  higher  energy  than  the  normal  radar  surveillance 
waveform.  This  principle  of  IRST  aided  phased-array 
radar  operation  leads  to  a  significant  increase  in  the 
confirmation  range  and  a  substantial  decrease  in  the 
radar  resources  required  for  track  maintenance.  As  a 
consequence,  more  tracks  can  be  maintained  and  more 
radar  resources  can  be  applied  to  track  initiation  [4, 
Sec.  14.9]. 

The  paper  explores  the  theoretical  bound  on  phased 
array  radar  allocation  when  tracking  highly  manoeu¬ 
vrable  anti-ship  missiles  (ASM)  using  a  collocated 
radar/IRST  sensor  combination.  The  assumption  is 
that  during  the  ASM  tracking,  the  phased  array  is 
requested  (to  allocate  a  beam  along  a  certain  direc¬ 
tion,  transmit  a  suitable  waveform,  and  process  the  re¬ 
ceived  echo)  whenever  the  track  error  exceeds  a  certain 
threshold.  The  motion  of  the  ASM  is  modelled  using 
the  quantized  acceleration  levels  in  three-dimensions. 
The  principal  tool  in  developing  the  bound  on  radar 
allocation  is  a  Cramer-Rao  type  error  bound  for  the 
estimation  of  linear  jump  Markov  dynamic  systems  [1]. 
Thus  given  a  dynamic  motion  model  of  an  ASM,  the 
IRST / radar  sensor  characteristics  and  a  tolerable  level 
of  target  state  estimation  error,  we  can  theoretically 
predict  the  maximum  average  update  time  required  for 
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the  phased-array  radar  (which  translates  into  the  re¬ 
quired  energy  per  unit  time).  The  significance  of  this  is 
that  one  can  predict,  even  before  the  system  is  built, 
the  upper  limit  on  the  capability  of  the  IRST/radar 
surveillance  system.  In  addition,  the  relative  mer¬ 
its  of  having  an  IRST  sensor  to  aid  in  tracking  can 
then  be  easily  quantified,  without  a  need  for  extensive 
Monte  Carlo  simulations.  The  model  of  the  phased  ar¬ 
ray  radar  that  we  consider  in  this  study  is  somewhat 
simplified  because  it  ignores  the  details  concerning  the 
radiated  waveform,  the  signal  processing  and  the  mul¬ 
tifunctional  capability  of  radar  (i.e.  the  interleaving 
of  functions  such  as  search,  track,  etc).  However,  the 
main  ideas  presented  here  can  be  directly  extended  to 
more  detailed  sensor  models. 

The  paper  is  organised  as  follows.  Section  2  presents 
a  mathematical  formulation  of  the  problem,  with  de¬ 
tails  of  the  ASM  motion  model  and  sensor  measure¬ 
ment  models.  Section  3  is  devoted  to  the  theoretical 
bound  for  radar  allocation.  Numerical  analysis  is  pre¬ 
sented  in  Section  4  and  the  conclusions  of  the  study 
are  drawn  in  Section  5. 


2  Problem  formulation 

We  consider  shipboard  collocated  phased-array  radar 
and  an  IRST,  targeted  by  a  highly  manoeuvrable 
ASM  (see  Fig.l).  The  two  sensors  are  assumed  to 
be  perfectly  registered.  The  target  (ASM)  state  is 
x  =  [x  x  y  y  z  z]T .  The  motion  model  that 
we  adopt  for  a  generic  ASM  is  described  next. 


Figure  1:  Illustration  of  the  engagement  scenario 


2.1  Dynamic  motion  models 

Various  target  dynamics  models  applicable  to  anti-ship 
missile  motion  (both  subsonic  and  supersonic)  have 
been  proposed  in  the  literature,  such  as  the  nearly  con¬ 
stant  velocity  model,  constant  acceleration  model,  co¬ 
ordinated  turn  model,  Singer  (coloured  process  noise) 
model,  “weave”  manoeuver  model,  “BUNT”  manoeu¬ 
vre  model,  to  name  a  few  [5,  6,  7].  Since  this  plethora 
of  (both  linear  and  non-linear)  models  is  neither  ex¬ 
clusive  nor  exhaustive,  we  adopt  an  alternative  mod¬ 
elling  approach  based  on  quantized  acceleration  levels 
[8,  9].  The  assumption  here  is  that  the  acceleration 
is  piece-wise  constant  and  its  mean  is  taken  from  a  fi¬ 
nite  set  of  acceleration  levels.  The  perturbations  upon 
this  constant  acceleration  is  modelled  as  white  Gaus¬ 
sian  noise.  Switching  from  one  acceleration  level  to 
another  is  modelled  by  a  first-order  Markov  chain. 

The  adopted  target  dynamic  model  is  expressed  as: 

xfc+i  =  Afcxfc  +  Bfea(rfe+i)  +  wfc  (1) 


where  k  is  the  index  assigned  to  the  continuous-time 
instant  tk, 
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Bfc 
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Here  I„  denotes  an  nxn  identity  matrix,  ®  is  the  Kro- 
neker  product,  Tk  =  tk+i  —  tk  is  the  sampling  interval; 
Wfc  is  a  6  x  1  vector  of  zero-mean  white  Gaussian  noise 
with  nonsingular  covariance  matrix 


Q  k=q-  block-diag  (0fc,  0fe,  0fc)  (4) 


where  q  is  a  parameter  related  to  process  noise  inten¬ 
sity  [7,  p.270],  and 


©fc 
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(5) 


The  acceleration  vector  a(rfc+i)  in  (1)  is  a  func¬ 
tion  of  a  discrete-valued  3D  random  vector  rfc+i  = 
[rx  ry  rz]T+1  which  determines  the  regime  (i.e.  ac¬ 
celeration  level)  of  target  motion  during  the  period 
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positive  integers.  The  accel- 

eration  vector  is  then 


a(rfc+i))  =  [a0rx  b0ry  c0rz]T  (6) 

where  do,  &o  and  Co  are  suitably  chosen  acceleration 
quantum  values  (or  units)  along  x,  y  and  z  axes,  re¬ 
spectively.  The  acceleration  is  thus  discretised  into 
S  =  (2 sx  +  1)(2 sy  +  1)(2 sz  +  1)  possible  levels.  The 
target  motion  regime  rk  can  switch  between  S  mod¬ 
els  in  a  random  manner.  The  evolution  of  the  motion 
regime  is  modelled  by  a  first-order  time-homogeneous 
Markov  chain  with  known: 

1.  transitional  probabilities 

TTij  =  P{rfc+i  =  j|rfc  =  i}.  (7) 

where  i  =  [ix  iy  iz]T,  j  =  \jx  jy  jz]T,  with 
i x '  jx  £  77..,. ,  iy ,  jy  £  TLy ,  i: z  •  j z  EE  A. 2 ,  the  t  ransi- 
tional  probability  matrix  [7Tij]  is  a  square  S  x  S 
matrix. 

2.  initial  regime  probabilities  pi(i)  =  P(ri  =  i), 
where  i  =  [ix  iy  iz]T ,  and  ix  £  1ZX,  iy  £  1ZV , 
i  z  £  7 lz. 

The  transitional  and  initial  regime  probabilities  are 
non-negative  and  normalised,  that  is: 

X^iCi)^1,  H7ri'j  =  1- 

i  j 


(8) 


2.2  Measurement  models 


The  IRST  provides  target  azimuth  and  elevation  mea¬ 
surements  at  regular  sampling  intervals  T.  The  loca¬ 
tion  of  the  IRST/radar  in  the  local  Cartesian  coordi¬ 
nates  is  given  by  vector  (xa.  yQ,  z0)T ,  however,  in  order 
to  simplify  notation  we  will  assume  that  xa  =  y0  = 
z0  =  0.  The  IRST  measurement  vector  at  time  k, 
zk  =  [0fc  efc]  is  then  given  by 

zfc  =  hfc(xfc)+vfc  (9) 


where 


hfc(Xfc)  = 


/  arctan  ^ 

xk 


arcsm 


y/xl+vl+zl 


(10) 


and  v[  is  the  IRST  measurement  noise,  assumed  to 
be  zero-mean  Gaussian  with  covariance  matrix  R/  = 
diag(ogj,  of  j). 

The  phased-array  radar  is  requested  to  provide  mea¬ 
surements  of  target  range,  range-rate,  azimuth  and  el¬ 
evation  based  on  the  adopted  radar  allocation  strat¬ 
egy.  The  radar  measurement  vector  at  time  k,  = 

[Pk  Pk  dk  efc  ] T  is  given  by 

Zfc=hfc(Xfc)+Vfc  (11) 


where 
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and  vj[)  is  the  radar  measurement  noise,  assumed  to 
be  zero-mean  Gaussian  with  covariance  matrix  = 
diag  {(T2p,(T%alR,alR). 

In  this  analysis  we  will  assume  that  the  radar  allo¬ 
cation  is  requested  whenever  the  target  track  error  (in 
position,  range,  elevation  or  azimuth)  exceed  a  prede¬ 
fined  threshold. 


3  Radar  allocation  bound 

The  key  idea  of  this  paper  is  to  compute  the  theoreti¬ 
cal  bound  on  required  radar  allocation  using  the  pos¬ 
terior  Cramer-Rao  lower  bound  (PCRLB)  for  target 
state  estimation.  Namely,  if  is  an  unbiased  target 
state  estimator  with  covariance  C  k,  then  the  following 
inequality  holds: 

Ck  =  E  [(xfc  -  xfc)(xfc  -  xfc)T]  >  Pfc  (13) 

where  P&  is  the  PCRLB  and  its  inverse  Jk  =  PjT  1  is 
the  Fisher  information  matrix  (FIM)  defined  in  [10]. 
Substantial  advances  have  recently  been  reported  in 
development  of  the  PCRLB  for  target  tracking,  such 
as  an  efficient  PCRLB  computation  in  the  context  of 
nonlinear  filtering  [11]  and  the  PCRLBs  when  measure¬ 
ments  are  of  uncertain  origin  [12,  13,  14,  15,  16].  In 
this  paper  we  will  use  the  best-fitted  Gaussian  (BFG) 
approximation  to  the  PCRLB  for  Markovian  switching 
systems  [1]  because  of  the  dynamic  model  of  the  ASM 
described  in  Sec. 2.1. 


3.1  The  BFG  approximation  for  jump 
Markov  linear  systems 

In  order  to  compute  the  PCRLB  for  the  dynamic  model 
described  in  Sec. 2.1,  we  have  to  augment  the  state  vec¬ 
tor  to  include  units  a o,  bo  and  cq .  We  point  out  that 
these  units  are  known,  hence  not  required  to  be  es¬ 
timated;  this  fact  will  be  reflected  in  the  initial  co- 
variance  matrix  and  the  process  noise  covariance.  We 
augment  the  state  vector  as  follows  [17]: 

■xk  =  [x  x  a0  y  y  b0  z  z  c0]T  (14) 


so  that  the  dynamic  equation  (1)  can  be  written  as: 


Xfe+i  =  Frkk+15ck  +wk 


where 


jprfc  +  l 
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Tk 

1 

0 


rj~\  2 

~2~ '  r ^,/c+i 
Tfc  *  i 

1 


(15) 


(16) 


Here  takes  values  x,  y  and  2  (in  this  order).  The 
process  noise  w/;:  is  again  zero-mean  Gaussian,  with 
covariance  matrix  Q&  =  q  ■  block-diag  ^©fc,  ©fc,  ©fc^, 
where  &k  =  diag(©fc,0). 

When  the  dynamic  equation  is  given  in  the  form  of 
(15),  we  can  apply  the  best-fitted  Gaussian  approxi¬ 
mation  to  compute  the  PCRLB  [1],  The  BFG-PCRLB 
computation  requires  us  to  specify  the  initial  target 
state  xD  and  its  covariance  matrix,  PQ,  the  initial 
regime  probabilities  Pi(i),  the  transitional  probabili¬ 
ties  7Ti j ,  process  noise  covariance  and  sensor  character¬ 
istics  (measurement  error  standard  deviations  and  the 
sampling  time). 

The  main  idea  of  the  BFG-PCRLB  approximation 
is  to  approximate  the  linear  jump  Markov  system  de¬ 
scribed  by  (15)  with  the  linear  non-switching  system 

xfc+i  ss  3>fcxfc  +  ufc  (17) 

where  u*,  is  an  “equivalent”  zero-mean  white  Gaus¬ 
sian  random  vector  with  a  covariance  matrix  X  k.  A 
straightforward  procedure  for  the  sequential  computa¬ 
tion  of  3>fc  and  Hk  ensures  that  the  first  and  the  second 
moment  of  models  (15)  and  (17)  are  identical  [1].  This 
procedure  is  given  in  Table  1. 

Example  1.  Consider  the  case  where  a  ship 
is  located  in  3D  space  at  (0, 0,  Ojm  and  the  ini¬ 
tial  location  of  the  ASM  is  (17321, 10000,  0)m.  The 
ASM  is  heading  towards  the  ship  with  the  speed 
of  400m/s.  The  initial  covariance  matrix  is  P0  = 
diag (aZ,a%,6,a%,aZ,6,a;,a%,5),  where  <rx  =  21  m, 
<jy  =  13.5  m  and  <rz  =  12.8  m  (computed  using  the 
first  (IRST,  radar)  associated  measurement  pair  and 
the  spherical-to-Cartesian  conversion).  Furthermore, 
av  =  10  m/s,  and  6  <C  1  is  a  very  small  value  which 
reflects  the  fact  that  the  quantum  values  ao,6o,co 
are  known  (d  has  to  be  non-zero  so  that  PQ  is  non¬ 
singular). 

The  transitional  probabilities  are  set  as  follows:  the 
non-switching  probability  is  7Ti  i  =  a  and  the  remaining 


Table  1:  The  BFG  computation  of  <1?^.  and  E*, 


Initialisation:  e  =  xD,  C  =  PD. 
For  k  =  1,2,3,... 


is  that  the  uncertainty  grows  with  the  larger  span  be¬ 
tween  the  minimum  and  maximum  acceleration  levels 
(in  this  case  red  corresponds  to  ±4 g  while  blue  corre¬ 
sponds  to  ±8<?) .  The  value  of  a  also  influences  the  size 
of  the  uncertainty  ellipsoid  (results  not  shown  here): 
the  higher  the  value  of  a,  the  smaller  the  uncertainty. 
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End 

1  —  a  is  equally  distributed  among  the  switching  cases 
so  that  normalisation  in  (8)  holds.  The  initial  non¬ 
accelerating  regime  probability  is  set  to  1,  i.e.  F(rXti  = 
0,ry,i  =  0,rZji  =  0)  =  1.  The  units  of  acceleration  are 
ao  =  bo  =  cq  =  1  g,  where  g  =  9.81m/s2.  Figure  2 
shows  the  2cr  uncertainty  ellipsoids  in  the  X  —  Y  plane 
of  the  BFG  approximation  for  a  =  0.9  and  q  =  0.001 
m2/s3  at  time  steps  tk  =  kTk ,  with  Tk  =  0.5s,  k  = 
0,1,....  In  red  we  show  the  results  for  sx  =  sy  =  4 
and  in  blue  for  sx  =  sy  =  8  (in  both  cases  sz  =0). 


3.2  FIM  measurement  contribution 

Once  we  determine  <3?/;  and  Efc,  we  can  apply  the 
Riccati-type  recursion  to  compute  the  FIM  [11]: 

Jfc  =  JP(k)  +  Jz(k)  (A:  =1,2,...)  (18) 

where  matrix 

Jp(fc)  =  (^-iJ^Li  +  Sfe-i)’1  (19) 

is  the  predicted  FIM  and  matrix  Jz(fc)  is  the  measure¬ 
ment  contribution  to  the  FIM. 

The  FIM  measurement  contribution  Jz(k)  in  gen¬ 
eral  depends  on  sensor  measurement  functions  (in  our 
case  h£  and  h^),  target  trajectory,  and  the  detection 
characteristics  of  sensors  (the  probability  of  detection 
Pd  and  the  probability  of  false  alarm  Pfa).  The  non¬ 
ideal  detection  characteristics  of  sensors  (Pd  <  1  and 
Pfa  >  0)  introduce  uncertainty  in  the  measurement 
origin,  which  reflects  itself  in  a  reduction  of  the  mea¬ 
surement  contribution  to  the  FIM.  A  comprehensive 
analysis  of  the  influence  of  Pd  <  1  and  P/a  >  0  on 
Jz(k)  is  presented  in  [16].  For  simplicity,  however,  in 
this  study  we  assume  zero  false  alarms  rates,  while  the 
effect  of  Pd  <  1  will  be  taken  into  account  as  follows: 


Figure  2:  Uncertainty  ellipsoids  of  the  BFG  ap¬ 
proximation  for  different  levels  of  quantized  ac¬ 
celeration.  In  red  is  the  case  TZX  =  lZy  = 
{-4ff,  -3g,  ...,0,  ...,3g,4g},  in  blue  is  Tlx  =  TZy  = 
{-8g,-7g,...,0,...,7g,8g}. 

The  first  observation  from  Figure  2  is  that  the  mean 
of  the  BFG  approximation  is  directed  towards  its  aim 
(the  ship),  despite  the  fact  that  the  ASM  is  manoeu¬ 
vring.  This  is  true  whenever  the  transitional  probabil¬ 
ity  matrix  is  symmetric  and  the  acceleration  space  is 
quantized  symmetrically  around  the  zero  acceleration 
in  all  directions.  The  second  more  obvious  observation 


J z(k)  =  Pd  ■  E  {H^R-'Hfc}  .  (20) 

The  resulting  bound,  referred  to  as  the  information  re¬ 
duction  factor  PCRLB  [15],  is  a  conservative  but  rea¬ 
sonable  approximation,  which  becomes  fairly  accurate 
after  a  few  initial  scans.  The  term  E  in  (20)  is  the 
expectation  operator  and  is  the  Jacobian  of  mea¬ 
surement  function  h^: 

Hfc(xfc)  =  [Vifch^(xfe)]T  (21) 

evaluated  at  the  true  value  of  the  x*;.  Depending  on 
the  source  of  a  measurement  at  time  k  (IRST  or  radar), 
function  hfc  takes  the  form  of  hjf.  or  h^,  respectively. 
Matrix  in  (20)  is  the  measurement  covariance,  and 
takes  the  form  of  R{.  or  R((,  depending  on  the  source 
of  the  measurement.  Similarly,  Pd  is  the  probability 
of  detection  of  a  sensor  (Pj  for  IRST;  Pj*  for  radar), 
and  typically  is  a  function  of  target  range.  If  at  time  fc, 
both  radar  and  IRST  measurements  are  requested,  due 
to  their  mutual  independence,  eq.(20)  takes  an  additive 
form: 

Jz(fc)  =  P(fE{(H«)T(Rf)-1H«} 

+  PiE{(H/fe)T(Ri)-1Hi}. 

Jacobian  is  a  4  x  9  matrix;  Jacobian  Hj[  has  the 
same  entries  as  the  lower  2x9  submatrix  of  H^.  Hence 
we  present  only  the  elements  of  H^,  which  are  easily 


obtained  by  differentiation.  The  non-zero  elements  of 
Hf  are  as  follows: 


Hf  [1, 1] 
Hf  [1, 4] 
Hf  [1, 7] 
Hf  [2, 1] 

Hf[2,4] 
Hfc  [2, 7] 

Hfc  [3, 1] 
Hf  [4, 1] 
Hf  [4, 4] 

Hf  [4, 1] 
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Pk 
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Xk 
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where  pk  =  \/x2k  +  yl  +  z\. 

In  order  to  compute  the  measurement  contribution 
to  the  FIM  as  in  (20),  we  would  need  to  average  the 
product  Hf  R/71Hfc  over  all  possible  realisations  of  the 
state  vector  xk.  In  general  this  can  be  done  numer¬ 
ically,  e.g.  via  a  sample  based  technique  [1].  In  our 
analysis,  however,  we  will  consider  one  particular  ASM 
trajectory  (see  Sec. 4.1)  for  which  we  will  compute  the 
measurement  contribution  to  the  FIM. 


3.3  Radar  allocation  requirement 

The  IRST  provides  its  measurements  at  a  regular  sam¬ 
pling  interval  Xj.  The  phased-array  radar  measurement 
is  required  for  track  update  whenever  the  predicted 
RMS  error  in  estimating  the  position  of  the  ASM  ex¬ 
ceeds  a  certain  threshold  77.  The  radar  allocation  test  is 
carried  out  every  Tk  <  Ti  seconds  and  formally  states: 
allocate  radar  if 

RMSE^f  «  >  r?, 

(22) 

where  Jp  is  the  predicted  FIM  of  (19).  The  exact  ex¬ 
pression  for  the  RMS  error  in  the  predicted  position 
would  involve  the  cross-terms  such  as  Jp  [1,4],  Jp  [1,7], 
etc.  (see  [10,  p.  9]  for  the  general  formulation  of  the 
PCRLB  for  the  nonlinear  transformation  of  the  state 
vector).  For  discussion  on  other  rules  for  radar  alloca¬ 
tion  see  [4,  Sec. 14. 4]. 


The  initial  speed  is  700  m/s,  with  acceleration  load  in 
the  horizontal  plane  oscillating  with  a  period  of  5  sec¬ 
onds.  The  height  of  the  missile  is  constant  at  10  m. 
The  trajectory  is  shown  in  Figure  3  (top-down  view) 
for  peak  acceleration  loads  of  4,  8, 12  and  16  g. 
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Figure  3:  Test  missile  trajectory  (with  varying  levels 
of  acceleration  loads) 

Numerical  analysis  is  carried  out  assuming  that  the 
radar  has  confirmed  the  initial  IRST  alert  and  that  the 
track  on  the  incoming  ASM  has  been  established  by 
fusing  the  initial  radar  and  IRST  measurements.  The 
elements  of  the  radar  measurement  covariance  matrix 
Rf  are  specified  as:  ap  =  30  m,  ay,  =  3  m/s,  ae,R  =  3.0 
mrad  and  =  4.5  mrad.  The  accuracy  of  IRST 
measurements  is  characterised  by  agj  =  0.25  mrad  and 
aej  =  0.64  mrad.  The  range  at  which  missile  tracking 
starts  is  16  km.  The  initial  covariance  matrix  is  for 
the  case  of  IRST/ radar  combination  specified  by  PQ  = 
diag(a-^ ,  a2 , 8,  cr2,  cr2 , 6,  a2,  a2, 5) ,  where  ax  =  21.8  m, 
cry  =  13.4  m  and  az  =  12.8  m,  av  =  10  m/s  and 
S  =  10 -8.  For  comparison,  we  also  consider  the  radar- 
only  (no  IRST)  case,  with  <jx  =  39  m,  ay  =  53  m  and 
az  =  75  m  (these  values  were  obtained  by  spherical- 
to-Cartesian  conversion) . 

The  probability  of  detection  Pd  for  both  sensors  is 
modelled  next.  In  the  case  of  the  IRST,  following  the 
description  in  [4,  Sec. 2. 3],  the  Pd  versus  range  curve  is 
shown  in  Fig. 4  in  red  solid  line.  This  curve  is  obtained 
assuming  Pd  =  0.9  at  the  nominal  range  of  17  km.  For 
the  case  of  a  radar,  an  attempt  has  been  made  to  incor¬ 
porate  the  effect  of  specular  multipath  on  Pd,  as  shown 
in  Fig. 4  in  blue  thin  line.  The  values  of  multipath  nulls 
for  the  radar  cannot  be  known  in  advance  as  they  de¬ 
pend  on  the  height  of  the  missile.  In  general  multipath 
also  deteriorates  the  accuracy  of  elevation  angle  radar 
measurements  [18],  although  this  effect  has  not  been 
included  in  the  model. 


4  Numerical  analysis 

4.1  Simulation  setup 

The  simulated  missile  trajectory  corresponds  to  a  su¬ 
personic  ASM  which  performs  a  weaving  manoeuvre. 


4.2  Numerical  results 

First  we  analyze  the  positional  RMS  error  as  a  function 
of  time  for:  (1)  the  case  where  an  IRST  is  available 
to  complement  the  phased-array  radar  in  tracking  an 
ASM;  the  IRST  is  scanning  the  horizon  with  T)  =  Is; 


Figure  4:  Probability  of  detection  versus  range  (Pj  in 
red;  Pj*  in  blue) 


(2)  the  radar-only  tracking  case.  The  results  are  shown 
in  Fig. 5  for  sx  =  sy  =  4,  sz  =  0,  a0  =  b0  =  c0  =  1  g, 
Tfc  =  0.2  s,  a  =  0.9,  q  =  0.0001  and  r)  =  100  m. 


Figure  5:  Positional  RMS  error  with  IRST  (blue)  and 
without  an  IRST  (red) 

From  Fig. 5  we  note  that  sharp  drops  in  the  RMS 
error  indicate  the  instances  when  a  radar  allocation 
was  requested.  Thus,  if  the  radar  is  operating  alone,  it 
is  requested  to  revisit  the  target  7  times;  with  addition 
of  the  IRST,  however,  the  radar  is  requested  only  3 
times. 

Next,  in  Fig. 6,  we  analyse  the  average  update  time 
of  the  phased-array  radar  as  a  function  of:  (a)  the  sam¬ 
pling  period  of  IRST  Tj  and  (2)  the  value  of  sx  =  sy 
which  corresponds  to  manoeuvre  ability  of  the  incom¬ 
ing  missile.  In  preparation  of  this  figure  we  used  the 
same  trajectory  and  parameters  as  indicated  above, 
unless  otherwise  stated. 

Fig. 6. (a)  was  obtained  for  sx  =  sy  =  8.  Observe 
that  the  IRST  sampling  intervals  (Tj)  of  0.5,  1  and 
2  seconds,  result  in  very  similar  radar  allocation  re¬ 
quirements.  There  is  a  dramatic  drop  in  average  Tr 
(i.e.  a  growing  demand  for  the  radar)  when  Tj  is  in¬ 
creased  from  2  seconds  to  8  seconds;  finally,  for  Tj  >  64 


(a) 


(b) 


Figure  6:  Average  radar  update  time  (in  seconds)  ver¬ 
sus:  (a)  the  IRST  sampling  period  Tj;  (b)  the  ma¬ 
noeuvre  ability  of  the  missile  in  the  horizontal  plane 
(sx  =  sy) 

seconds,  the  IRST  becomes  irrelevant  and  the  system 
performs  as  if  the  radar  is  alone. 

Fig.6.(b)  was  obtained  for  sz  =  0  and  Tj  =  Is.  Ob¬ 
serve  that  the  higher  the  manoeuvre  ability  of  the  mis¬ 
sile,  the  smaller  the  average  update  time  of  the  radar 
(i.e.  the  greater  demand  for  the  radar). 

4.3  Comparison  with  an  EKF 

Next  we  compare  the  theoretical  bound  on  average 
radar  sampling  interval  Tr,  with  the  update  time  of 
an  Extended  Kalman  filter  (EKF),  designed  to  process 
both  radar  and  IRST  measurements  using  the  appro¬ 
priate  measurement  functions  and  hj[.  The  EKF 
assumed  a  constant  velocity  dynamic  motion  model, 
with  a  large  amount  of  white  process  noise  to  account 
for  ASM  manoeuvres.  The  radar  allocation  test  is  car¬ 
ried  out  as  in  (22),  except  that  J~ 1  is  replaced  by  the 
EKF  predicted  covariance.  The  radar  update  times  of 
the  EKF  were  averaged  over  10  Monte  Carlo  runs  and 
plotted  against  the  IRST  sampling  period  in  Fig. 7. 

Over  the  range  of  IRST  sampling  periods  examined, 


Figure  7 :  Average  radar  update  time  (in  seconds)  ver¬ 
sus  the  IRST  sampling  period  fy:  blue  line  is  the  the¬ 
oretical  (upper)  bound,  green  line  is  the  EKF  result 

the  EKF  achieved  only  about  half  of  the  upper  bound 
on  average  radar  update  time.  It  may  be  possible  to 
approach  the  bound  using  a  better  tracking  filter,  such 
as  an  IMM  [7,  9]  with  a  number  of  acceleration  states. 

Note  that  the  results  do  not  take  into  account  the 
potential  errors  in  associating  the  measurements  to  the 
track.  If  measurements  were  incorrectly  associated,  a 
decrease  in  the  average  radar  update  time,  or  equiva¬ 
lently  an  increase  in  radar  resources,  would  be  required 
to  maintain  the  desired  track  error  limit. 

5  Conclusions 

The  paper  presented  a  theoretical  analysis  of  phased- 
array  radar  allocation  requirements  for  tracking  anti¬ 
ship  missiles  with  an  integrated  radar/IRST  surveil¬ 
lance  system.  The  main  tool  in  the  analysis  was 
the  posterior  Cramer-Rao  lower  error  bound  for  tar¬ 
get  tracking,  and  therefore  the  obtained  results  for 
the  average  radar  update  time  indicate  only  the  up¬ 
per  bound.  The  presented  analysis  allows  us  to  quan¬ 
tify  the  IRST  benefits  in  the  anti-ship  missile  defence, 
without  a  need  for  extensive  Monte  Carlo  simulations. 
The  analysis  was  carried  out  using  simplified  radar  and 
IRST  models,  characterised  by  measurement  accura¬ 
cies  and  probabilities  of  detection.  Future  work  will 
extend  the  analysis  by  (1)  using  more  realistic  sensor 
models  and  (2)  calculating  a  more  accurate  PCRLB 
in  cluttered  environments,  via  the  enumeration  of  the 
measurement  sequences  as  in  [16].  Finally  it  will  be 
necessary  to  verify  the  theoretical  results  by  compari¬ 
son  with  the  experimental  data. 
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